Modified equation of state, scalar field, and bulk viscosity in Friedmann universe 
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A generalized dynamical equation for the scale factor of the universe is proposed to describe the 
cosmological evolution, of which the ACDM model is a special case. It also provides a general 
example to show the equivalence of the modified equation of state (EOS) and a scalar field model. 
In the mathematical aspect, the EOS, the scalar field potential V(ip), and the scale factor a(t) 
all have possessed analytical solutions. Such features are due to a simple form invariance of the 
equation inherited which determines the Hubble parameter. From the physical point of view, this 
dynamical equation can be regarded as the ACDM model with bulk viscosity, an existence content 
in the universe. We employ the SNe la data with the parameter A measured from the SDSS data 
and the shift parameter 1Z measured from WMAP data to constrain the parameters in our model. 
The result is that the contribution of the bulk viscosity, accumulated as an effective dark energy 
responsible for the current cosmic accelerating expansion, is made approximately ten percent to that 
of the cosmological constant. 

PACS numbers: 98.80.Cq, 98.80.-k 
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I. INTRODUCTION 



The cosmological observations have provided increas- 
ingly convincing evidence that our universe is undergoing 
a late-time cosmic acceleration expansion |l|, |2| ■ In order 
to explain the acceleration expansion, physicists have in- 
troduced a new fluid, which possesses a negative enough 
pressure, called dark energy. According to the observa- 
tional evidence, especially from the Type la Supernovae 
0, and WMAP satellite missions we live in a fa- 
vored spatially flat universe consisting approximately of 
30% dark matter and 70% dark energy. The simplest 
candidate for dark energy is the cosmological constant, 
but it has got the (in-)famous and serious fine-tuning 
problem, while the also elusive dark matter candidate 
might be a lightest and neutral supersymmetry particle 
with only gravity interaction. Recently, a great variety of 
models are proposed to describe the universe, partly such 
as scalar fields, exotic equation of state, modified grav- 
ity, inhomogeneous cosmology model and super horizon 
backreaction 6j . It is also instructive to investigate by 
concrete models the possibility that the both mysteri- 
ous dark components, dark matter and dark energy, may 
be two facets of a single complicated fluid, that is, they 
originate from a unified dark fluid in the universe. From 
Einstein equation we know that the right hand side is 
intimately related to the geometry or gravity, and on the 
other side the energy-momentum tensor also involves a 
quantum field vacuum. According to Ref . , it is tempt- 
ing to investigate the properties of cosmological models 
starting from the equation of state (hereafter EOS) of 
dark energy directly and by testing whether a given EOS 
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is able to give rise to cosmological models reproducing 
the available dataset. We hope the situation will be im- 
proved with the new generation of more precise observa- 
tional data upcoming. 

The observational constraints indicate that the cur- 
rent EOS parameter w = p/p is around —1 3], quite 
probably below —1, which is called the phantom region 
and even more mysterious in the cosmological evolution 
processes. In the standard model of cosmology, if the 
w < — 1, the universe shows to possess the future finite 
singularity called Big Rip (3, E} ■ Several ideas are pro- 
posed to prevent the big rip singularity , like by introduc- 
ing quantum effects terms in the action [l(J,[ll[. Based on 
the motivations of time-dependent viscosity and modified 
gravity, the Hubble parameter dependent EOS is consid- 
ered in Ref. [3, [l2( ■ The Hubble parameter dependent 
term in this EOS can drive t he p hantom barrier being 
crossed in an easier way Il3| . Different choices of 
the parameters may lead to several fates to the cosmo- 
logical evolution [14J . Recently, the equivalence between 
the modified EOS, the scalar field model, and the modi- 
fied gravity is demonstrated in Ref. 0], with a general 
method to calculate the potential of the scalar field for a 
given EOS presented. 

In Refs. 0, 0, 0] , the bulk viscosity in cosmology 
has been studied in various aspects. Dissipative processes 
are thought to be present in any realistic theory of the 
evolution of the universe. In the early universe, the ther- 
modynamics is far from equilibrium, the viscosity should 
be concerned in the studies of the cosmological evolu- 
tion. Even in the later cosmic evolution stage, for exam- 
ple, the temperature for the intergalactic medium (IGM) , 
the baryonic gas, generally is about 10 4 K to 10 6 K and 
the complicated IGM is rather non-trivial. The sound 
speed c s in the baryonic gas is only a few km/s to a 
few tens km/s and the Jeans length A yields a term as 
an effective viscosity c s A. On the other hand, the bulk 
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veloci ty o f the baryonic gas is of the order of hundreds 
km/s 19]. So it is helpful to consider the viscosity ele- 
ment in the later cosmic evolution. It is well known that 
in the framework of Fredmann- Robertson- Walker (FRW) 
metric, the shear viscosity |20j has no contribution in the 
energy momentum tensor, and the bulk viscosity behaves 
like an effective pressure. Moreover, the cosmic viscos- 
ity here can also be regarded as an effective quantity as 
caused mainly by the non-perfect cosmic contents inter- 
actions and may play a role as a dark energy candidate 

m 

In this letter, we show that in the framework of Fried- 
mann universe, the general EOS 

P = (7 - !)/> + Po + w hH + w H2 H 2 + w dH H, 

corresponds to a scalar field model whose potential has 
got the form 

V(ip) = V (e-^ + Cie- pv/2 + C 2 ). 

We will present analytically that the equation for the 
scale factor derived from the above EOS or scalar field 
model is more general than the ACDM model can show, 
and it has possessed an exact solution. The natural inter- 
pretation of this model is involved to the bulk viscosity. 
Concerning on the different forms of the bulk viscosity co- 
efficient, we propose three parameterized H(z) relations 
and use the observational data to constrain the parame- 
ters. 

This paper is organized as follows: In the next section 
we give the generalized dynamical equation for the scale 
factor and show a transformation to reduce the dynami- 
cal equation of the scale factor a(t) into a linear differen- 
tial equation. In Sec. Ill we demonstrate that the EOS is 
corresponding to a scalar field whose potential can be ex- 
actly solved. In Sec. IV we find that there exists a form 
invariance related to the variable cosmological constant 
which satisfies a renormalization equation. In Sec. V we 
use the SNe la data with the parameters A from large 
scale survey and shift 1Z from cosmic background radia- 
tion data to constrain our model. Finally, we present our 
conclusions and discussions in the last section. 



II. VISCOUS DARK FLUID DESCRIBED BY 
AN EFFECTIVE EOS 

We consider the FRW metric in the flat space geometry 
(fe=0) as the case favored by WMAP data 

ds 2 = -dt 2 + a(t) 2 {dr 2 + r 2 dfl 2 ), (1) 

and assume that the cosmic fluid possesses a bulk viscos- 
ity £. The energy-momentum tensor can be written as 



T lu , = pU li U v + (p-S0)H ia ,, 



(2) 



where in comoving coordinates = (1,0,0,0), and 
Hfj, u = 9fj,u + U^Uu 22]. By defining the effective pres- 
sure as ft = p — t^9 and from the Einstein equation 



Rfj,v — \911vR = STrGTpv, we obtain the Friedmann equa- 
tions 



8ttG 
A-kG 



(p + 3p). 



(3a) 
(3b) 



The conservation equation for energy, T.°", yields 

p+( P + ft)8 = 0, (4) 

where 9 = U!i = 3d/ a. 

In our previous work [l4j . we have considered the fol- 
lowing EOS form with the same notations 

p = (7- l)p + p a + w H H + w H2 H 2 + w dH H, (5) 

where H is the Hubble parameter and w x are their cor- 
responding coefficients. We have obtained the exact so- 
lution to the scale factor and showed that 
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C = Co + Ci- + C 2 t 
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is equivalent to the form by using the above EOS. By 
defining 
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the dynamical equation of the scale factor a(t) can be 
written as 



a 
a 



37 - 2 a 2 



1 d 1 



(11) 



The analytical solution for 7 ^ is given out as 
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For the case 7 = 0, the solution is 



a(t) = oq exp 



1 



)(< 



(t-io)/Ti 



(12) 



Tift -t )' 
T 2 



(13) 

The five parameters in Eq. (JSJ) are condensed to three 
free parameters in Eq. (flip or its solution, later the best 
fit analyses in Sec. V enable us to obtain the physical 
evolution of the universe. 
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It is interesting that there exists a transformation 

y = a 3 ^ /2 (14) 
to reduce Eq. to a linear differential equation of y(t) 



37 _ 



1 



(15) 



which can be solved easily. The equation of a(t) can be 
written as a more general form 



37 -2 a 2 la 1 
2 a? + J\a + T, 



C 



,37/2 • 



(16) 



and we can also use the transformation of Eq. (|14f) to 
obtain the linear differential equation of the function y(t) 



The interesting case is that u(<j>) and V((j>) are deter- 
mined in terms of a single function f(<j>) as 

^) = 4/'W. y(^) = 4(3/W 2 + /'W)- (23) 

One can check that the special solution cj> = t, H = f(t) 
satisfies the scalar-field equation. The following relations 
is obtained in order to solve the function f(cj>) for a given 
EOS, 



By defining a new field ip as 



<P = / d<t>y/\w(4)\, 



(24) 



(25) 
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(17) 



the action can be rewritten as 



We can check directly that the following transformation 



y = a ^/ 2 + CT% 



(18) 



can reduce the nonlinear Eq. (|16[) to its corresponding 
linear Eq. JT^. If 7 = 4/3, i.e., w = 1/3, Eq. flty 
becomes 



a a 2 1 a 1 C 
a a 2 T\ a a? ' 



(19) 



which can be interpreted as the radiation dominated uni- 
verse when the curvature of the universe is concerned. 



III. SCALAR FIELD AND MODIFIED 
GRAVITY 

Starting from the action for the gravitational with the 
matter fields, we can show that an EOS for the universe 
contents corresponds to a scalar field model. Generally, 
for a given EOS, the potential V(<p) often has got no 
analytical solution. Here we demonstrate that the corre- 
spondent scalar field model for the EOS of Eq. JSJ luckily 
has possessed an exact solution. First, we revisit the gen- 
eral procedure proposed in Ref. ^5| to relate a scalar field 
model for a given EOS. 

Starting from the action of the scalar-tensor theory 



5 / d^xj—g (^R - 1^)0^4, - V(</>) 



the energy density and the pressure are 



(20) 
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P = —LJ 
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2 +V(0), P =^LU 



Vicji). (21) 



Combining the above equations and the Friedmann equa- 
tions, one obtains 



W (0)</> 2 



rH, v{4>) 



1 



(3H 2 + H). 



(22) 



S = 



jRT^d^-V^)). (26) 



The energy density and the pressure is now given by 

p = ±^ 2 + V&), P = ±^ 2 -%). (27) 

Now we summarize the procedure. For a given EOS, 
by using Eqs. 124(1 . one can obtain an equation for f(4>); 
then solving the equation gives /(</>). Using Eqs. I|23|l . 
one obtains lo(4>) and V(<p). And by employing Eq. I|25|) 
to transform the variable <fi to if, finally, one obtains the 
V((f). In Ref. [3, several examples are presented such as 
p = wp, and we show in this section that a more general 
form of Eq. © can also have had an analytical solution 
of the potential. 

By using the Friedmann equations, Eq. (JSJ) can be 
rewritten as 



P= (7-l)p- 



The corresponding equation for f((f>) is 



f{4>) = f/W>) 2 - 
The solution of this equation is 

f{<j)) = acoth(-^-a<p 

where 



« 2 T 2 ' 
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(29) 



(30) 



(31) 



Then we obtain the uj(4>) and V((j)). The integration of 
Eq. I|25|) gives the relation between cj> and ip as 



<P = — V3|tI In 

K 



tanh(37a0/2) 



tanh(37a0o/2) 



(32) 
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FIG. 1: The w-t relation of the quintessence case (panel A) 
and the phantom case (panel B). 



which gives 

coth(ya^) = coth(-^-"<? ! 'o)exp(--V3|^|¥'). (33) 

So substituting Eq. l|3Tl|) to Eq. l(2*3*)l and using the above 
equation to transform the variable <j> to ip, we obtain 



V(<p) = ^ coth" ( ^-acPo 



3(2-7) 



.JL e xp(-|,/»> 



1 



+ 3 A 

3*/ 2 T?a 2 2 



34) 



This potential has got the form 

V(<p) = V Q (e-^ + de-^' 2 + C 2 ) . (35) 

As a special case, p = (7 — l)p is given out in Ref. [l5j . 
Taking the limit of Ti — > 00 and a — > of Eq. and 
omitting the constant term, we obtain 



where 



V(cp) = V Q e- K V ItIv, 
2(2-7) 



K 



3 7 ^ 2 ^- 



(36) 



(37) 



To see the dynamical effects of the scalar field, by defin- 
ing the EOS parameter w = (±p 2 - V)/(±<p 2 + V), the 
evolution of w with respect to the time t is illustrated in 

Fig.cn 

The action for a modified gravity in Einstein frame 



5' 



d X\/—g ( TT~2 R + ^matter 



(38) 



is shown to be related to a modified EOS [T5|. In 
Ref. 0, the following two equations are derived in the 
framework of modified gravity, 
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= =H 2 
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,> - f(R) + 6 ( H + H 2 - H-) f'(R), 



(39) 



= —{2H + 3H 



2 —H - 3H Z 



J 1 P + f(R) 

d 2 d 



dt 2 



2H-)f'(R) 7 (40) 



where R = 6H + 12H 2 . If we adopt the approximation 
H « H 2 [12, and assume f(R) = foV~R, by combining 
Eqs. and we can obtain 

p=( 1 -l)p-^f 0l H-\ 1 H 2 : (41) 

which is a special case of Eq. (JjJJ . Additionally, the scalar- 
tensor theory in Einstein frame is mathematically equiv- 
alent to the modified gravity in Jordan frame 15]. The 



makes the ki- 



conformal transformation g^ v 
netic term in the action of scalar-tensor theory vanish, 
then one obtains the Jordan frame action 



S = 



d A x, 



-9 



2k 2 



-R 



<p)\. (42) 



If the equation of motion of ip can be solved as p = 
<p(R), one obtains the action of the modified gravity in 
Jordan frame. Thus, the Hilbert-Einstein action with an 
additional term f(R) has made effective contributions 
similar to those caused by a scalar field, as well as a 
modified EOS. 



IV. VARIABLE COSMOLOGICAL CONSTANT 
AND RENORMALIZATION GROUP EQUATION 

The model of variable cosmological constant is another 
alternative to explain the cosmic evolution, in order to 
overcome the serious fine-tuning problem. In the frame- 
work of variable cosmological constant model, the evolu- 
tion of the scale factor is determined by both the Fried- 
mann equations and the renormalization group equation 
(RGE) for the cosmological constant [24l |25| , 



37 - 2 d 2 37 + I 



A, 



<2A 



-ap 2 M 2 



(43) 
(44) 
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Here we have already used the EOS p = (7 — l)p to 
eliminate the p and p in Friedmann equations. With 
the choice of the renormalization scale p = H |24|, the 
variable A is determined by 

' /A 1 (45) 



d\nH 4tt 2 



aH 2 M 2 . 



The solution is 



A(i) = A + t[H(t) 2 -H(t ) 2 ]M 2 , 



(46) 



which has got the form of A = Co + C2H 2 . Compared 
Eq. (g3J) with Eq. JTTJ), we find that if Eq. (gSJ is formly 
invariant under the transformation 



p, — > p + Sp, 



(47) 
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the solution of RGE becomes 



k = C + C 1 H + C 2 H 2 . 



(48) 



Substitute this result to Eq. I)43|) . we obtain a equation 
which has the form of Eq. 1)11(1 . Especially, if 5/i = 0, 
then C\ — 0. It is very interesting that concerning on the 
bulk viscosity, modified EOS, scalar field model , modi- 
fied gravity, and the variable cosmological constant can 
be described in one generally dynamical equation which 
determines the scale factor. 

Actually, the evolution equation of the Hubble param- 
eter, 



rp2 ' 
J 2 



(49) 



has possessed a form invariance for H — > H + Hq, i.e 
37 



1 



H =-T H ' + {T 1 - 2Ho ) H+ n 



i 



(50) 



It is this form invariance that gives several interesting 
features of the model, such as both the a(t) and V(ip) 
have analytical solutions for the general EOS of Eq. J5J. 



V. A NEW CONTENT OF THE UNIVERSE 
AND DATA FITTING 



In Eq. (HU 

mological constant 



l/Tf plays the role of the effective cos- 



If To 



oo, the H — z relation is 



H(z) = H [n m {i + zf' 2 + (i - n m )] 



(51) 



We proposed a parameterized H-z relation as the follow- 
ing form 

H 2 = H 2 [Vl m {\ + z)^ + Vl v (l + zfl 2 + \-n m -VL v l (52) 

where Sl v = 2/{Z ! YTiH ) and l-Q m -fi„ = 2/(3^H^). 
Note that Eq. (| 1 1 ft can have several interpretations |14| . 
such as the " inflessence" model [2(j when the first term 
and third term alternatively dominate. The three terms 
in the right hand side of Eq. (|llfl are proportional to 
H 2 , H 1 , and H°, respectively. In the early times, the 
first term is dominant, which may lead to inflation if 
7^0. In the medium times, the second term domi- 
nates, which leads to deceleration if T\ < 0. In the late 
times as current stage, the third term is dominant, which 
leads to cosmic acceleration expansion behaving as the de 
Sitter universe if Tf > 0. So, a single equation may de- 
scribe three epoches of the cosmological evolution. In 
this paper, however, we only focus on the interpretation 
of bulk viscosity for this model, as the viscous universe 
has been discussed for various cosmology evolution stages 
with very naturally physical motivations. And the vis- 
cosity and the dissipative processes in describing physical 
universe have been studied in various aspects, for exam- 
ple nn m mom. 



TABLE I: 


w 


universe content 


contribution to H\z) 


—1 


vacuum 


(l + zf 


-1/2 


(effective) viscosity 


(l + 2 ) 3/2 


-1/3 


curvature 


(1+^) 2 





dust 


(1+^) 3 


1/3 


radiation 


(1 + *) 4 


1 


stiff matter 


(1 + ^) 6 



So far, the universe contents and their dynamical con- 
tributions are listed in Table |U Compared with the 
ACDM model, we use three cosmological scenarios as pa- 
rameterizations of the H-z relations, which are listed in 
Table EH 



• (i) The ACDM model, the simplest model to ex- 
plain the dark energy. 

• (ii) The viscosity model without cosmological con- 
stant. 

• (iii) The bulk viscosity is constant, so that the bulk 
viscosity has got the dynamical effects of (1 + z) 3 / 2 . 

• (iv) The bulk viscosity has the form of £ = £o + 
(jd/a, where the constant term has the dynamical 
effects of (1 + z) 3 / 2 , and the term proportional to 
H has an effect to change 7 to 7. (See Eq. (7J and 
Ref. [3) 

The observations of the SNe la have provided the di- 
rect evidence for the cosmic accelerating expansion for 
our current universe. Any model attempting to explain 
the acceleration mechanism should be consistent with the 
SNe la data implying results, as a basic requirement. The 
method of the data fitting is illustrated in Refs. |32l l33| . 
The observations of supernovae measure essentially the 
apparent magnitude m, which is related to the luminosity 
distance by 



m(z) = M + 5 log 10 D L (z), 



(53) 



where Dl{z) = (Hq/ 'c)cIl(z) is the dimensionless lumi- 
nosity distance and 



Also, 



d L (z) = 



M = M 



c(l + z) 
Ha 



1 



E(z' 



■dz'. 



Bl ° V IMpc J 



25, 



(54) 



(55) 



where M is the absolute magnitude which is believed to 
be constant for all supernovae of type la. We use the 157 
golden sample of supernovae la data compiled by Riess et 
al. to fit our model. The data points in these samples 
are given in terms of the distance modulus 



fiobs{z) = m(z) - M obs (z). 



(56) 
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FIG. 2: The la (solid line), 2<r (dashed line), and 3a (dotted 
line) contour plots of model (iii). 



FIG. 3: The la (solid line), 2a (dashed line), and 3a (dotted 
line) contour plots of fl m -Q v relation in model (iv). 



The x 2 is calculated from 



E 



Hobs(zi) - M! - 51og 10 D Lth (zi;c a ) 



1 2 



&obs(Zi) 



A-0A69\ 2 . fTZ-1.716\ 2 



0.017 



0.062 



I ■ (57) 



where A4' = M. — M b s is a free parameter and 
DLth(zi',c a ) is the theoretical prediction for the dimen- 
sionless luminosity distance of a SNe la at a particular 
distance, for a given model with parameters c a . The pa- 
rameter A is defined as yjl 



A = 



0.35 



0.35 



0.35 



dz 



£(0.35) J E(z) 



1/3 



and the shift parameter 1Z is [35, 36] 

rZ " dz 



K 



'fir. 



15(z) 
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FIG. 4: The la (solid line), 2a (dashed line), and 3a (dotted 
line) contour plots of Q. m -w relation in model (iv). 



VI. CONCLUSION AND DISCUSSION 



We will consider the ACDM model for comparison and 
perform a best-fit analysis with the minimization of the 
X 2 , with respect to M' , fi m , fl v , and 7. The results are 
summarized in Table ITT1 From the results, we see that 
the bulk viscosity part has made approximately 10% con- 
tributions to that of the cosmological constant. In model 

(ii) , the cosmic acceleration expansion is due to the bulk 
viscosity, without the cosmological constant, and with 
comparing this model is rather disfavored. For model 

(iii) , Fig. [3 plots the likelihood contour of the param- 
eters fl m and fl v . We can see that the bulk viscosity 
contributes approximately 10% of the cosmological con- 
stant. For model (iv), Fig.[3|and Fig.|Ijplot the likelihood 
contours oftt m — n v and Q m ~w (w = 7— 1), respectively. 



In conclusion, we have shown that several different ap- 
proaches to explain the current accelerating universe ex- 
pansion can give the same form of a dynamical equation 
for a(t). Also in the sense that different terms in the 
right hand side of the dynamical equation may domi- 
nate correspondingly different periods we call the media 
described by the general EOS a Dark Fluid. The case 
of the ACDM model can be regarded as a special case 
of this equation descriptions. In the framework of Fried- 
mann universe, the correspondences of the bulk viscosity, 
EOS, scalar field, modified gravity, and variable cosmo- 
logical constant are summarized in Table ITTT1 The scalar 
field model as a prototype is intensely used to study the 
dynamical behaviors for the scalar factor in the litera- 
ture because of its simplicity in formulation and treat- 
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TABLE II: Summary of the best fit parameters 





Models 


best fit of parameters 


x 2 


(i) 


H 2 = Ho[Q m (l + a) 3 + 1 — f2 m ] 


Q m = 0.283 


177.84 


(") 


h = Hi{n m (i + z)V 2 + i-n m ] 


fl m = 0.435 


319.78 


(iii) 


h 2 = Hi[n m (i + z f + + z) 3/2 + 1 - n m - n„] 


(n m ,0,) = (0.281, 0.065) 


177.38 


(iv) 


# 2 = #S[n ro (i + ^ + ^) 3/2 + 1 - fim - a,] 


(O ro , n„, 7) = (0.298, 0.053, 1.004) 


177.32 



merit; here we have provided a general and exact solu- 
tion of the scalar field which can be either quintessence 
or phantom. Mathematically, this equation enjoys sev- 
eral interesting features, such as the EOS, the potential 
V((p) of the scalar field model, and the scale factor a(t) 
all have got analytical solutions. This equation can be 
transformed to a linear equation and has possessed a form 
invariance related to the Hubble parameter. Physically, 
the bulk viscosity is a most natural interpretation of this 
model. We have proposed three cosmological scenarios 
which contain the contributions from the bulk viscosity 
and used the 157 golden sample of SNe la data with the 
parameter A and TZ to constrain our model. The spe- 
cial model that only the bulk viscosity contributes the 
cosmic accelerating expansion without the cosmological 
constant may be excluded. And the (effective) bulk vis- 
cosity has made approximately 10% contributions to that 
of the cosmological constant. 

The puzzling cosmic dark components: dark matter 
and dark energy, responsible mainly for large scale struc- 
ture formation of universe and cosmic accelerating ex- 
pansion as well as our universe evolution fate as we now 
understand in the standard hot big bang and inflation 
models have challenged our previous intelligence on the 
physical world. A unification picture description for the 
two elusive dark composites either from complicated fluid 
dynamics, modified gravity, inhomogeneous cosmology 
or quantum field models with introducing more degrees 
of freedom and supported by more precious experiments 
like LHC and PLANCK in 2007 is certainly valuable for 
us to uncover the mysterious dark side of the Universe, 
which even will bring us new knowledge on fundamental 
physics. 



Note: In this brief presentation we find by numerical 
calculations that parameters A and TZ respectively from 
large scale survey and cosmic background radiation de- 
tections significantly affect the fitting results for some 
models. For example, if we do not use A and TZ, the 
minimum \ 2 of Model (ii) is 179.36, which is acceptable. 
Furthermore, the best-fit results of Model (iii) are listed 
in Table llVl From this table, it is obvious that the re- 
sult by using the 157 golden samples (hereafter 157 as in 
the below table) with parameters A and 1Z is very dif- 
ferent from that by only using the 157 golden samples, 
for which we are trying to figure out the physical reasons 
in late study. It is very likely that we can obtain more 
accuracy results with high level certainties by next year 
(2007) PLANCK mission, continuous large scale struc- 
ture and later large sample of SNe la survey, to which we 
are very confident. 



TABLE IV: Fitting results for the model H 2 = # 2 [fi m (l + 
z) 3 + Q V (1 + z) z/2 + 1 - Q m - Q v ] with different data 









x 2 


157 


0.330 


-0.583 


174.53 


157 + A + 71 


0.281 


0.065 


177.38 


157 + .4 


0.300 


0.063 


177.36 


157 + ft 


0.303 


0.045 


177.34 


A + TZ 


0.316 


0.117 


1.88 x 10" 


A 


0.283 


0.031 


1.34 x 10" 


K 


0.266 


0.033 


3.47 x 10" 
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